Cooling electrons by magnetic-field tuning of Andreev reflection 
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A solid-state cooling principle based on magnetic-field-driven tunable suppression of Andreev 
reflection in superconductor/two-dimensional electron gas nanostructures is proposed. This cool- 
ing mechanism can lead to very large heat fluxes per channel up to 10^ times greater than cur- 
rently achieved with superconducting tunnel junctions. This efficacy and its availability in a two- 
dimensional electron system make this method of particular relevance for the implementation of 
quantum nanostructures operating at cryogenic temperatures. 
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In recent times, the vast interest in mesoscopic sys- 
tems brought thermoelectric phenomena under the spot- 
light [ij in the context of solid-state refrigeration 
Promising solid-state refrigeration schemes operating at 
sub-Kelvin temperatures rely on superconducting mi- 
crostructures [3|] and were recently shown to yield re- 
markable electron [3] and phonon [^| temperature reduc- 
tion. In such systems, a normal metal (N) region is cou- 
pled to a superconductor (S) through an tunnel insulat- 
ing barrier. Quasiparticle cooling occurs because of the 
existence of the superconducting gap, which allows only 
the more energetic electrons to escape from N. In this 
scheme the characteristics of the insulating barrier are 
crucial: on one hand it must suppress Andreev reflec- 
tion (in order to limit NS current), on the other it 
must be sufficiently transmissive to maximize hot-carrier 
tunneling into S. Despite the remarkableperformance of- 
fered by this refrigeration mechanism [J, |3, 0, 0| , insu- 
lating barriers suffer from sub-optimal transmissivity at 
energies above the superconducting gap thereby limiting 
the efficiency of heat extraction. 

Recently some of us Q proposed to suppress An- 
dreev reflection while keeping highly-transmissive barri- 
ers by employing ferromagnets rather than insulators. In 
this Letter we present a different operational principle 
and investigate heat transport in superconductor/two- 
dimcnsional electron gas (2DEG) nanostructures in the 
presence of a localized magnetic field. The high trans- 
missivity available with this configuration yields a heat 
current per channel that can exceed by up to a factor 
of 10^ what is currently achievable with insulating bar- 
riers. This can lead to efficient electron refrigeration, 
even in macroscopic 2DEG regions. The origin of this 
effect stems from the fact that a localized magnetic field 
applied next to a S/2DEG interface can strongly sup- 
press Andreev reflection while affecting only marginally 
the above-gap transmissivity [§]. Furthermore, the fine 
modulation of the localized field, together with the tun- 
ability of the electron-phonon interaction in 2DEGs, al- 
low full control of the heat current for the optimization 
of this refrigeration principle. 
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FIG. 1: (color) (a) Scheme of the proposed setup. Ferromag- 
netic strips (F) placed on top of the structure generate at each 
S /2DEG junction localized magnetic fields whose strength can 
be controlled by changing the orientation of magnetization 
(At), (b) Idealized semiclassical picture of the Andreev re- 
flection suppression mechanism: an electron (e) impinges on 
the localized magnetic field (yellow region) at the critical an- 
gle (Scrit) above which no Andreev reflection takes place. 



The system we propose consists of a 2DEG region in- 
serted between two identical S electrodes and is sketched 
in Fig. [Ija). Ferromagnetic strips (F) are placed on 
top of each S/2DEG junction and a voltage V is applied 
across the structure. Junctions are located at a; = x;, 
where i = 1,2 corresponds to the left and right contact, 
respectively, d long, t thick F layers are positioned sym- 
metrically on top of each junction, zq above the 2DEG 
[see Fig. [IJa)]. F layers are assumed to have homoge- 
neous magnetization Al. Heat and electron transport 
in the 2DEG region are affected by the perpendicular 
magnetic field produced by the x-component of the mag- 
netization, M.^ — |Al|cos0. The 0-component of the 
magnetic field (B) can be controlled by rotating Al by 
means of a weak magnetic field ^ . For the sake of clarity, 
we shall first focus on a single S/2DEG ballistic contact, 
namely the left one: we set xi = and z = in the 
middle of the ferromagnetic layer. When t <^ zq and 
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t d, the z-component of the magnetic field generated 
by the left F strip in the 2DEG plane is given by [I^, 11 1 
S,{x) = B^[JC{x+!i,zo)~ICix-lzo)]9ix), where 
IC{x,z) = —zd{x'^ + z'^)^^, /io is the vacuum magnetic 
permeability, and 0{x) is the Heaviside function [l^ . 
This expression for Bz {x) assumes that the magnetic field 
is entirely screened in the superconductor. In reality, the 
magnetic-field profile generated by the F strip shows a 
second peak located deep in the S region, which could 
possibly suppress superconductivity locally, but which 
does not affect Andreev reflection at the S/2DEG in- 
terface Q . The mechanism leading to Andreev reflection 
suppression relies on the effect of the localized magnetic 
field on the orbital motion of the quasiparicles. A sim- 
ple semiclassical picture, although idealized, captures the 
essential features of our system. As shown in Fig. [Ijb), 
the Lorentz force acting on electrons (e) impinging on 
the superconductor and on Andreev-reflected holes (h) 
introduces a critical angle (^crit) above which incoming 
electrons cannot be Andreev reflected. 

Within the Landauer-Biittiker scattering approach, 
the heat current J{Vc) through the junction biased at 
voltage Vc can be written as 

1 f°° 

J=T dsY^i^- eV,) {[AA-(e) - 7^^(e)] /(e - eV,) 

-7^X (£)/(£ + eV,) - [M^e) - TZ'^ie) - 7^1(e)] /(e)} ,(1) 

where TV*^ is the number of open channels for spin-(7, 7?.^ 
{TZ"^) is the spin-dependent Andreev (normal) reflection 
probability, and / is the Fermi function. The scatter- 
ing probabilities can be evaluated numerically through 
a recursive Green's function technique based on a tight- 
binding version ^ISJ of the Bogoliubov-de Gennes equa- 
tion. In the absence of spin-flip scattering the latter reads 



Wix) A{x) 
A*(.t) -H"{x) 



^l — <y 



(2) 



withH- = {v-eK{x))^^{v-eA{x)) + V,{x)+U{x)- 

, where p = —ih V, A(x) = Ay{x)y is the vector 
potential in the London gauge, and u {v) is the coher- 
ence factor for electron- (hole)-like excitations of energy 
e, measured from the S chemical potential e|. The po- 
tential tl{x) describes subband mismatch between S and 
the 2DEG. The Zeeman splitting in the semiconductor 
is given by Vij{x) = ^crg* iibB{x), where /iB is the Bohr 
magneton, and g* is the 2DEG effective g-factor. The 
quasiparticle mass and the S pairing potential are given, 
respectively, by m(x) = m*9{x) + mo9{—x) and A(x) = 
A9{—x), where nie {in*) is the free-electron (effective 
2DEG) mass. Within the tight-binding description, H'^ 
and A are matrices with elements {Ti.)ij = oJiSij — "yS^ijy 
and {A)ij = AoSij , where uji is the on-site energy at site i, 
and the 7 is the hopping potential ({•••} stands for first- 
nearest- neighbor sites) |1J|. The presence of the mag- 
netic field is included in the hopping potentials through 
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FIG. 2: (color) (a) Heat current J of a ballistic S/2DEG 
system extracted from the 2DEG side vs bias voltage Vc ap- 
plied across the junction, calculated for several angles (j) at 
T = 0.3 Ao/fce. /cb is the Boltzmann constant. Dash-dotted 
line represents J for a N/2DEG junction at <j!> = 0°. (b) J 
vs Vc calculated for several temperatures at </> = 0°. (c) J 
calculated at the optimal bias voltage across the junction vs 
T at = 0°. J is maximized at T = Topt ~ 0.25 Ao/fcB- (d) 
Efficiency {rj) of the cooling process calculated at the opti- 
mal bias voltage vs T at = 0°. The junction normal state 
resistance is i?N(<^ = 0°) = 1.07 x lO" Q.. 



Peierls phase factors [15|, thus taking into account the 
possibility of Landau level quantization in the 2DEG. 

Figure 2(a) displays how a strong enhancement of heat 
transfer off the 2DEG region occurs by decreasing (j) from 
(j) = 90°, i.e. by increasing the perpendicular magnetic 
field intensity. This is due to the increasing suppression 
of Andreev reflection 6| at the junction induced by the 
magnetic fringe field [8| [see also Fig. [Hb)]. Notably, 
J is positive (implying heat extraction from the 2DEG) 
at (f) = 60°, and is maximized around (j) — 40° where it 
reaches the value J~9xl0~^ Ag/e^i?N corresponding 
to - 0.5 pW/um-i, with R^ic/) = 0°) 1.07 x lO^Sl the 
junction normal state resistance. The smooth and grad- 
ual variation of J for 0° < ^ < 60°, demonstrates that a 
precise tunability of the heat flux is possible. At each an- 
gle, J is maximized by an optimal Vc- For 0° < (p < 40° 
this optimal bias is ~ 0.8 Ao/e. The behavior of J vs 
for several temperatures at ^ = 0° is shown in Fig. 2(b). 
At each temperature, J is maximized by an optimal bias 
voltage which, at the lowest temperatures, is 14 ~ Ao/e. 
Figure 2(c) shows J vs temperature, calculated at each 
optimal bias voltage for = 0°. In particular, J is 



maximized ai T = T, 



opt 



0.25Ao/fcB where it reaches 
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FIG. 3: (color) (a) Sketch of the system considered: the 2DEG 
is assumed to be in quasiequilibrium at temperature Tc, and to 
interact with the lattice phonons at temperature Tl through 
piezoelectric coupling (Pl)- Thermal load due to Joule dissi- 
pation (Pj) in the 2DEG, as well as the heat flux (J) flowing 
through the system are indicated with additional arrows, (b) 
System equivalent electric circuit: the 2DEG is shown as a 
lumped resistor (Pioad), and V is the total bias voltage ap- 
plied across the symmetric structure. 



Jmax — 6 X 10^^ Ap/e^iiN- It is interesting to compare 
this S /2DEG cooling scheme with the more conventional 
metal/insulator/superconductor (NIS) junction configu- 
ration. For the S/2DEG system, the junction normal- 
state transmissivity is T — 2iTh{e'^ R^My)^^ — 4.4 x 
10~^, where A/p = 54 is the total number of open chan- 
nels at the Fermi energy, and the maximum heat current 
is J,„ax ~ 9 X 10-2 A2/e2i?i^ [see Fig. 2(a)]. Consid- 
ering a state-of-the-art AI2O3 tunnel junction with spe- 
cific resistance Rs = lOOfi/Ltm^ l^, i5j, correspon ding to a 
normal-state transmissivity T^unnei — 5 x 10~^ [1^, the 
maximum heat flux is (Jmax)Nis — 6 x IS.'^/e^R^ 
at T — O.BAo/Zeb [ill- Comparing the two cases, one 
can readily obtain that ~ ^ T^^ 1 (^^7^) — 

1.2 X lO** ^ "^p" ^ i-e., the cooling power per channel in 

the S/2DEG is up to 4 orders of magnitude larger than 
currently attainable with insulating barriers. The differ- 
ence in the heat-transport properties originates entirely 
from the magnetic- field control on the S/2DEG contact 
which, in contrast to an insulating barrier, strongly sup- 
presses the sub-gap conductance while only marginally 
affecting the junction normal-state resistance. 

An important figure of merit of any refrigeration pro- 
cess is the efficiency r]{Vc) = J(yc)[VcIiVc)]~^ [H, where 
I is the electric current flowing through the junction. The 
efficiency rj versus T calculated at each optimized bias 
voltage for = 0° is displayed in Fig. 2(d). Efficiency 
values as high as 25% can be obtained at T ~ 0.25 Ao/fce. 



These results can be readily transferred to the case 
here of interest, electronic refrigeration in S/2DEG/S 
structures biased at voltage F [see Fig. 1(a)]. We focus 
on the quasiequilibrium limit [l|, where strong electron- 
electron interaction thermalizes the 2DEG at an effec- 
tive electron temperature To that may differ significantly 
from that of the lattice Tl [see the scheme in Fig. 3(a)]. 
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FIG. 4: (color) (a) Electron temperature Te in the 2DEG of 
a S/2DEG/S refrigerator vs V calculated for several angles (p 
at Tl — 0.3 K. (b) Tc vs V for different lattice temperatures 
Tl at </> = 0°. (c) Te vs V for several charge densities of the 
2DEG (n2DEG) at = 0° and Tl = 0.3 K. 



This limit can be attained by assuming a large enough 
I [see Fig. Ufa)] so that 1^ < I <^ II, with 1^ {II) 
the electron-electron (electron-phonon) scattering length 
19|. At low lattice temperatures (typically below 1 K) 
20|, piezoelectric electron-phonon scattering is the main 
mechanism that transfers energy into III-V 2DEG sys- 
tems fill 22 , [2^ , and this regime will be considered in 



the following. We model for simplicity the semiconductor 
region as a lumped element of resistance i?ioad connected 
to the junctions [see the system equivalent circuit in Fig. 
3(b)]. The stationary Tg reached by the 2DEG is ob- 
tained by solving 



2J(14, Te, Tl) + PL(Te, Tl) - Pj(K, Te, Tl) 
2Vc + Rio^dliVc, Te, Tl) - = 0, 







(3) 



where the first equation requires that the sum of all ther- 
mal fluxes in the 2DEG is zero, while the second is Kirch- 
hoff's voltage law for the whole system [see Fig. 3(b)]. 
Moreover, PL(Te,TL) is the energy-exchange rate due to 
piezoelectric coupling (24*], Pj = Pioad-f (K;, Te, Tl)^ is the 
Joule dissipation in the 2DEG, and the factor 2 in the 
first of ^ accounts for the fact that there are two iden- 
tical junctions. We implicitly assumed that the electron 
temperature of the S electrodes equals Tl [H. 

In the following [see Fig. 1(a)] we choose parameters 
typical of Ino.75Gao.25As 2DEGs [25*]. Figure 4(a) shows 
the electron temperature Te in the 2DEG vs bias voltage 
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V as obtained from the solution of ^ , for several angles 
(/) at Tl = 0.3 K. By reducing cj), can be drastically low- 
ered with respect to Tl. The refrigeration process turns 
out to be finely tunable for < (f) < 50° with a minimum 
Te ~ 50 mK at y =i 2Ao/e for </> < 40°. The elec- 
tronic temperature Te for different lattice temperatures 
at = 0° is displayed in Fig. 4(b) (Tl values coincide 
with the To value at ^ = 0). This cooling principle makes 
it possible to reach a minimum Tc around 12 mK starting 
from Tl = 0.2 K (temperature reduction 95%). Even 
for lattice temperatures as large as Tl = 0.5 K, Te re- 
ductions of 42% can be achieved. As can be inferred 
from the piezoelectric energy-exchange rate Pl(To,Tl) 
(see Ref. ^24]), .?^piczo(T) can be changed by varying the 
charge density in the 2DEG (due to the fcp ^-dependence), 

— 1/2 

thus leading to a ng ' -dependence of Pl. This is a rele- 
vant peculiarity of 2DEGs as compared to normal metals. 
Figure 4(c) shows Te vs V at </> = 0° and Tl = 0.3 K cal- 
culated for several values of the charge density (n2DEG) 
in the range ~ 1.1 x 10^^ — 1.75 x 10^^ m~^. Changes in 
the minimum Te of about 30 mK (at V ~ 2 Aq / e) can be 
obtained by varying n2DEG- 

In summary, we have analyzed heat transport in 
S/2DEG nanostructures, showing that a large heat flux 
per channel can be achieved in the presence of localized 
magnetic fields. The precise tunability offered by mag- 
netic (with fields as low as some mT [9|) as well as elec- 
trostatic control gives enhanced freedom for the control 
of electronic temperatures. Materials as luxGai-ajAs (for 
x > 0.75) [13] or InAs ^ 2DEGs (providing Schottky 
barrier-free contacts with a metal) in combination with 
Al and Co appear as ideal candidates for applications 
ranging from fully-tunable refrigerators to sensitive radi- 
ation detectors [l| operating at cryogenic temperatures. 
Partial financial support from MIUR under FIRB pro- 
gram RBNEOIFSWY is acknowledged. 
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